New renormalization procedure for perturbed ordinally differential equations with 
given initial condition and prescribed secular terms 
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For perturbed ordinally differential equations, a procedure of renormalization group method is 
proposed. To uniquely obtain renormalized solutions for given initial conditions, the procedure 
assumes that the extra integral constants yielded by perturbative expansion are determined as they 
cancel out all the prescribed secular terms. A sufficient condition for the renormalized solutions 
coinciding with exact solutions is also stated. 
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To obtain solutions to perturbed differential or dis- 
crete equations, a perturbative expansion is usually ap- 
plied. But naive perturbation often yields secular terms 
due to resonance phenomena, and the secular terms pre- 
vent us from getting approximate but global solutions. 
The renormalization group method § is one of singu- 
lar perturbation techniques, and gives the approximate 
but global solutions to a wide variety of differential [|[ ||] 
and discrete |{| equations. The method removes the 
secular terms by renormalizing integration constants con- 
tained in a solution to an unperturbed equation. Tempo- 
ral evolution of the renormalized constants is governed by 
a renormalization group equation, which is mathemati- 
cally regarded as an envelope equation for a family of 
naive solutions |J. 

However, a perturbative solution up to the fc-th order 
of a small parameter e contains (k+ 1) times as many in- 
tegration constants as the exact solution does because of 
the presence of undetermined constants appearing from 
the homogeneous solution from each order, and hence we 
cannot obtain the unique renormalized solution from a 
given initial condition without any assumptions. In this 
letter, we introduce an assumption, and propose a sys- 
tematic procedure to uniquely obtain the renormalized 
solutions both for oscillating naive solutions and for non- 
oscillating ones. We study perturbation up to the first 
order of e throughout this letter. 

Consider the perturbed differential equation of the 
form 



dx 

n = m 



■sg(x), x e K'\ 



(1) 



where x is a n-dimensional column vector, e is the small 
parameter, i.e. |£|<cl, and the unperturbed part of this 
equation, dx/dt — f(x), is integrable. To obtain a naive 
solution, x is expanded as a positive power series of e; 



,(o) 



■ ex 



(i) 



The equation of the zeroth order, 



the unperturbed equation, is 



dt 



= /(* (0) ), 



a;( ) = (f)(W(t;A), A £ R n , where each element of the n- 
dimensional column vector A is an integration constant. 
We assume the function <jy^ is smooth with respect to 
both t and A. Denoting the Jacobian of / by Df, the 
first order equation is expressed as 



dxW 
dt 



Df(xM)xM +g(x^). 



(2) 



The general solution to Eq.(g|) is composed of the ho- 
mogeneous solution and a special solution. The homoge- 
neous solution has extra n integration constants, which 
are written as the n-dimensional column vector a, and is 
denoted by (jP-> (t; A, a) accordingly. The homogeneous 
solution can be rewritten as 



(pW(t;A,a) 



dA 



(t; A)a, 



(3) 



where dcjr- Q > /dA is the Jacobi matrix. To show this equa- 
tion, we note that the function (f>^ (t; A + ea) also sat- 
isfies the first order equation: {d<jy- Q > / dt){t\ A + ea) — 
f((f>(°\t;A + ea)). Comparing the left- and the right- 
hand-sides in terms of the first order of the Taylor se- 
ries, we find that {d(t> {0) /dA) a is the general solution 
to the homogeneous part of Eq.(||). The special solu- 
tion is the sum of an arbitrary homogeneous solution, 
cf>^(t; A, a'), and the inhomogeneous solution, which has 
no extra constants and is denoted by ip(t;A). Since 
the homogeneous solution is not affected by the per- 
turbation g(x), secular terms are included by ip(t;A) if 
they exist. We therefore decompose it into two ports: 
ip(t;A) = i/>W(t;A) + ^ ns \t- A), where ^ s \t;A) and 
(/i' ns '(t;i) consist of secular and nonsecular terms, re- 
spectively. In the end, the naive but general solution up 
to the first order is formally written as 



„NS 



and the general solution to this equation is denoted by 



(t) = 0f°3 (t; A) + e[0 {1) (t; A, a) + (t; A) + (t; A)} , 

(4) 



where we have replaced a + a' with a. 

The first and crucial step for constructing the unique 
renormalized solution is to determine a local solution 



2 



from an initial condition, x NS (t ) = Xo, where the lo- 
cal solution must be valid around to. Let A and a be 
the values of A and a determined by the initial condition, 
respectively. The number of elements of A Q and a is 2n, 
while the initial condition can determines values of n el- 
ements among them. We thus introduce an assumption 
that the extra integration constants, a, cancel out all the 
secular terms of the naive solution (0) at t = to'- 



<j ) < - 1 \to;Ao 1 ao)+iP is \to;Ao) = 0. 



(5) 



This condition uniquely determines the value ao from Aq 
and to unless accidental degeneracy occurs. The local 
solution being valid around to, then, becomes 



x NS (t;t ) = ^°\t;A ) 



e[^ {1 \t-A ,ao) 
+ ^ s Ht;A )+^ s \t;A )}, 



(6) 



and the value Ao is uniquely determined from the initial 
condition x NS (t ;t ) — xo with Eq.(||). 

The second step is to connect the local solution 
x NS (t; to) with another local solution being valid around 
/I. The time [i may be arbitrarily chosen from a suitable 
neighborhood of to where the local solution x N (t;to) is 
valid, and hence we regard the integration constant A as 
a function of fi, and denote it by A(/j,). Similarly, the 
extra integration constant a is regarded as a function of 
A(/i) and \x, and is denoted by a{A{p),\x). We remark 
that A(to) — Aq and a(A(t ),to) — ao- Then the local 
solution around /i is similar to Eq.(|^): 

x NS (t; /i) = (t; A(»)) + efrW (t; A(fi),a(A(n), M )) 
+ ^ s \t;A(v))+^ ns Ht;A(n))}, 

and a(A(fj,), ji) must be satisfied the similar condition to 
Eq.©: 

0« (/i; a(AOi), M )) + ^ (s) (m; A(/i)) = 0. (7) 

Using Eq.(|^), the values of the extra integration con- 
stants are formally determined as 



This is the renormalization group equation, and must be 
solved with the initial condition A(to) — Ao- We remark 
that the renormalization group equation is determined by 
the extra integration constants through the dependence 
on A and /i, which is assumed as Eq.(ph. 

Finally, the renormalized solution, x^ (t), is obtained 
by evaluating the naive solution x NS (t;to) at to = t and 
A(to) = A(t), because we are interested in the solution 
being valid around the time t. The renormalized solution 
is, therefore, expressed as 



x^{t) =x Nb {t-t) = 4>^{t;A{t))+s^ ns \t;A{t)), 

(9) 

by using the condition (^). 

Let us apply this procedure to the following three ex- 
amples. We use the same symbols introduced above, and 
denote elements of x, A and a by Xj, Aj and aj, j = 1, 2, 
respectively. 

Example 1 - Harmonic oscillator 

The first example is a harmonic oscillator with linear 
perturbation: 



_d fxi 
dt \x 2 



1 

-1 





-xi 



First, we diagonalize the matrix found in the right-hand- 
side of this equation with the diagonalizing matrix R and 
the coordinate transformation which are 



R 



1 i 

i 1 



R 



respectively. The original equation is transformed as 



±fyi 

dt U 2 



i 
-i 



£ I iyi- 2/2 
-yi - m 



(10) 



The naive solutions up to the first order are calculated 

as 



a(A(^) j/i )=-(^0x;A( / i))) V (s W(m))- (8) 



The connection condition, 



(*;<o) 



satisfied up to the first order of e by putting A(fi) 
A(to) + 8 A and deciding the 8 A as 



8 A 



-e[a(A(to),n)-a(A(to),t )}. 



Remembering 8 A = A(/j,) — A(to) and taking the limit 
/i — > to, we obtain the relation between differential co- 
efficients at to- dA/d[i(to) = —eda/d[i(A(to),to)- Since 
the initial time to may be chosen arbitrarily, we regard 
the relation as the differential equation of the form 



dA 
dfi 



s^(A(ti),n). 



y» s (t) = A 1 e it + e 
y? s (t) = A 2 e- it + 



ai e lt + -AxtJ 1 + \-.A 2 e- lt 
2 4i 



a 2 e-n - -A 2 te- U ~ —A^ 1 



4i 



(11) 



and <j>^ and are determined as 



4> {1) {t-A,a) = 



a 2 e~ 



^\t;A) = - 



A-yte* 
-A 2 te~ u I 



respectively. From the condition (Q), the extra integra- 
tion constants, a% and a 2l are determined as 

ai(A(p),fi) = -iAt(fj,)fj,/2, a 2 (A{fi),n) = iA 2 (n){j,/2, 

(12) 
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and the renormalization group equations are 

dA 2 



dAi i 



d/i 



Solving these renormalization group equations with the 
initial condition (Ax(tp), A 2 (to)) = (Axfi,A 2< o) , which 
is determined by Eqs.(|ll|),(12) and the initial condition 
(l/i(*o)j Ste(*o)) = (1/1,0,2/2,0), the renormalized solutions 
are uniquely obtained as 



Ax e fe(*"*o)/2 e it + e -fa(t-to)/2 e -« 

4i 



4/ 



^i,o e 



ie(*-i )/2. 



Note that the renormalized solutions satisfy the original 
equation (|l^) up to the first order, and the frequency 
of the renormalized solution, 1 +e/2, approximates the 
frequency of the exact solution, i/l + e, up to the first 
order of the Taylor series accordingly. 

Example 2 - Jordan cell 

The second example is an equation having a Jordan cell 
in the unperturbed part: 



dt \x 2 



1 






-xx 



(13) 



The naive solutions up to the first order are calculated 
as 



„NS 



x 2 



(*) 
(*) 



Ax 
A 2 



Ant 



a 2 



ax 
Axt 



a 2 t 
1 



A,C 



The solutions of the zeroth order do not oscillate, and 
no resonance occurs accordingly. Consequently, it is not 
obvious which are secular terms, but we regard the last 
two terms in each of Xi S {t) and x 2 s {t) as the secular 
terms since their absolute values exceed absolute values of 
zeroth order when et > 1. To vanish all the secular terms 
simultaneously at t — fj,, the extra integration constants, 
ax and a 2 , must be determined as 



ax(A(fj,),(J,) 
a 2 (A(n),fi) 



-A!0V/2-A 2 0V/3, 
A 1 (ji)+A 2 (ji)ii 2 /2. 



Then the renormalization group equations are 



— = £(fj,Ax+fi J A 2 ), 



dA2 (A 



fiA 2 



We emphasize that the renormalized solutions, x^{t) = 
Ax(t) + A 2 (t)t and x 2 G (t) = A 2 (t), satisfy the original 
equation ( [l3"| ) with the aid of the renormalization group 
equations, and hence they are exact solutions. 

We remark that an equation which has a Jordan cell in 
its unperturbed part is also discussed in Ref.|7[, but one 



of the renormalized integral constants is selected from a 
homogeneous solution of the first order. 

Example 3 - Nonlinear equation 

The third example is an equation whose unperturbed 
part is nonlinear: 



d_ fx x 
dt \x 2 



3'2 



-3xj/2 



if 

2 l-exp(-2xi] 



(14) 



This is Einstein's equation for a Friedman-Robertson- 
Walker universe with dust, and x\(t) is the logarithm 
of the scale factor of the universe and e is the sign of 
the spatial curvature. The naive solutions up to the first 
order are 

x» s {t)=\n{t + Ax) 2 ^+A 2 

ax(t + Ai)- 1 + a 2 - jj-(* + Ai) 2 / a exp(-2A 2 ) 



„NS 



(t) = -(t + AxY 



-ax(t + AxY 
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[t + Ai)- 1 ' 3 exp(-2A 2 ) 



The solutions of the zeroth order do not oscillate again, 
and hence secular terms are not determined automati- 
cally. If we decide the secular term as 



^ S \t-A) 



-9{t + Ax) 2 / 3 exp(-2A 2 )/20 




then the renormalized solution coincides with one ob- 
tained in Ref . B with changing the independent variable 
as t = (fi + Ai) 2 / 3 . The renormalized solution a;^ G (t) 
improves the naive solution Xx S {t) and reproduces ex- 
panding and contracting feature of the exact solution. 
On the other hand, if we decide the secular terms as 



^ s \t;A) 



-9(t + Ax) 2 ^ 3 exp(-2A 2 )/20 \ 
-3(t + Ax)~ 1/3 exp(-2A 2 )/10J ' 



then the renormalization group equations are 



cL4i 

d/i 
dA 2 

d/i 



e±e- 2A >^ + Ax) 2 / 3 , 



- £ 2 6 



-2A- 



: (M + ^i)- 1/3 , 



and the renormalized solutions satisfy the original equa- 
tion (|l4|) exactly. 

In Examples 2 and 3, renormalized solutions become 
exact when nonsecular terms do not appear in the first 
order. In the following Theorem and Corollary, we show 
that the feature holds in general. As a preparation of the 
Theorem, we define the functions 7 = s, ns, as 



gW(<t>W(t;A)) 



<9V (7) 
dt 



{t-A)-Df(4> {0) {t-A))^\t-A), 



4 



and the symbol <SC/(x )s>2 as the higher orders of the 
Taylor series of f(x RG ), i.e. 



«/(l KC )»2 = 



k=2 



where represents the A:-th order derivative of /. 

Theorem: Let x RG (t) be a renormalized solution to the 
ordinally differential equation ([j]) up to the first order. 
Then the following relation holds: 



dx 



RG 



dt 



(t)-[f(x RG (t))+sg(x RG (t))] 



OA y ' w 7 V dA 
(cf>W ( t; A (t))) - £> 9 (0(°) (t; A(i)))V> ( " s) (t; A(t)) 
-«/(a; i?G )»2 - e«:9(x flG )» 2 



(15) 



Proof: In this proof, we omit arguments of functions 
which are evaluated at fi = t and A = A(t). Note 
that, from the definitions, the relations hold: 84>^ /dt = 
f((f> {0) ), dip (ns) /dt = D/(^°))Vj ( ™ s) + .g ( " s) O (0) ) and so 
on. 

Derivating the renormalized solution (^|) with respect 
to t, we obtain the following equation 



dx RG 



00<°> da 



d/ JXT ' "OA dfi 

D/(^ (0) )V (ns) +g (ns) (0 (O) ) 



94 9/i 

Next, partially derivating Eq.(^) with respect to we 
find the equation of the form 

where we used Eqs.(^) and (0). From the two equations 
mentioned above, we obtain the following relation 



dx RG 
dt 



fd^ r 



\ OA 



\ dA ) 



On the other hand, from Eq.^), the Taylor series of 
f(x RG (t)) +eg(x RG (t)) becomes 

f(x RG ) + eg{x RG ) = /(0 (o) ) + eg{^ 0) ) + eDf{^)^ 
+ e 2 Ag(0 (O) )V' ( " s) +«/(^ G )»2 + e« 5 (x flG )» 2 . 

Consequently, we obtain the relation stated in the theo- 
rem. 

Corollary: The renormalized solution up to the first or- 
der is the exact solution to the original equation if the 



naive solution of the first order does not include nonsec- 
ular terms. 

Proof: From the assumption, ^( ns ) = 0, and <<C 
f{x RG {t)yz>2 =*mg{x RG {t))^>2 = accordingly. Conse- 
quently, the right-hand-side of Eq.(|l^) vanishes, and the 
renormalized solution x R (t) exactly satisfies the original 
equation. 

As summary, we propose a procedure for uniquely con- 
structing a renormalized solution from a given initial con- 
dition. The basic idea of the procedure is to determine 
extra integral constants as they cancel out all the secular 
terms simultaneously at the initial time. As a result of 
the procedure, the theorem is stated that a renormalized 
solution is the exact solution to an original equation if 
the naive solution does not include nonsecular terms in 
the first order of a small parameter. Note that renormal- 
ization group equations are not always autonomous even 
if original equations are. 

The renormalization group method is also available to 
obtain reduced equations as the renormalization group 
equations ^ [l(], pl| . The procedure proposed in this 
letter clarifies the relation between the renormalization 
group equation with the way how we decide secular 
terms, and hence it seems useful not only obtaining ap- 
proximate solutions but also reducing original equations. 
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